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The realization of phase discontinuities across metasurfaces has led to a new class of reflection and 
refraction. Here we present theory and experiment on the discontinuous propagation of wavepackets across 
subwavelength-thickness meta-atoms. Using acoustic waves, we observe the process of wavepackets 
traversing a meta-atom with abrupt displacements, which appear as path discontinuities on a space-time 
diagram. We construct a tunable meta-atom from two coupled resonators at —500 Hz, map the 
spatiotemporal trajectories of individual sonic pulses, and reveal discontinuities at the meta-atom where the 
pulses exit at a time —50 ms ahead or behind their arrivals. Applications include thin acoustic metasurface 
lenses. 



Metamaterials are artificial structures engineered on a sub -wavelength level to produce exotic wave 
propagation characteristics 1 " 5 . Ranges of constitutive parameters exhibited by their bulk forms have 
been extended to beyond classically available ranges in both electromagnetism and acoustics: single and 
double negativity as well as a wide range of double-positive parameters have brought many remarkable phe- 
nomena into reality, including negative refraction, reverse Doppler effect, superlensing, and cloaking 6 " 9 . 

Recently a new direction for metamaterials research has emerged by the introduction of metasurfaces 10-25 . 
Holloway et at. 10 define a metasurface as any array of metamaterial unit-cells (meta-atoms) on a surface for which 
the thickness and periodicity are small compared to the wavelength in the surrounding media. For the simplest 
case of one-dimensional (ID) waves propagating in a waveguide, a metasurface is represented by a dot, or single 
meta-atom 10 . With significantly smaller loss, metasurfaces produce abrupt phase and amplitude changes, enab- 
ling new mechanisms for refraction, polarization conversion and beam shaping 15 " 21 . 

Despite intensive investigations on abrupt phase shifts, few studies have been reported on the explicit action of 
metasurfaces on wavepackets, there having been only one reported case of the direct observation of wavepackets 
passing through a metasurface: Dolling et al. 26 measured an abrupt advancement of the group time-shift across an 
electromagnetic metasurface in the near-infrared. However, a detailed picture of the transmission process in the 
vicinity of the metasurface was not obtained. Here we present the first spatiotemporal measurement of a 
wavepacket as it propagates across a meta-atom. To achieve this we use an acoustic system consisting of a tunable 
meta-atom made up of two coupled resonators. 

Geometry and theory of the acoustic meta-atom. We start by briefly reviewing pulse transmission through 
a dispersive medium or an impedance boundary with dispersive transmission characteristics 27 , showing 
analytically that metasurfaces producing a phase discontinuity for monochromatic waves also engender an 
abrupt displacement for passing wavepackets. This analysis applies to both electromagnetic and acoustic 
waves. The phase shift 3 across a metasurface exhibits a strong dispersion. Therefore, the beat pattern 
resulting from the generic case of the superposition of two monochromatic waves of nearby frequency is 
altered by the difference in phase shift as they traverse the metasurface. Figure 1 shows schematically the case 
of a higher frequency wave at angular frequency co h encountering greater phase advancement 3 than a lower 
frequency wave at a>i [see Fig. la]. In this case it is straightforward to show that the temporal advance of the 
envelope x g [see Fig. lb] is given by x g = A3/A(o. 
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Figure 1 | Schematic diagram illustrating the origin of an abrupt 
discontinuity at a metasurface by means of the simultaneous 
transmission of two monochromatic waves of different frequencies a?i 

and c&fr (a) Here, the higher frequency component pi experiences more 
phase advancement than the lower frequency component p 2 as the waves 
pass through the metasurface, resulting in (b) a temporal shift in the beat 
pattern. The blue dotted curve in (b) represents two frequencies with the 
same phase shift. The red solid curve in (b) represents the case when the 
higher frequency component experiences more phase advancement. 

For the more general case of a wavepacket of central frequency oo c 
and amplitude pi (t) = m (t) cos (co c t) = m(t)Re{exp( — ico c t)} at the 
position of the metasurface, where m(t) is a slowly varying envelope 
function, the complex form of the transmitted wavepacket p T (f) is 



GO' 

p T (t) = | P I (w)T(w)e- i0Jt dw, 



where Pj (co) is the Fourier transform of pi (t) and T(co) is the 
complex transmission coefficient. If the spectral width of the wave- 
packet Aco is sufficiently narrow, T(co) can be approximated to first 
order as T(w) ~ |T(co c )| exp(h p co c — h g (w — co c )), where i p = 3(co c )/ 
oo c and x g = ddldoo are respectively the phase shift and the group 
time-shift evaluated at the frequency oo c . Here we are assuming that 
the transmission amplitude does not vary much over the spectral 
width of the pulse, Aco. Substituting for T(co) in the above equation, 
we have 



~p T {t) = \T(w c ) 



-z'co(t+t. 



g) dco. 



The transmitted wavepacket is therefore given by 

p T (t)=Re{p T (t)} = \T(oj c )\m(t + T g )cos(oj c (t + T p )). 



Clearly, the envelope of the exiting wavepacket has the same shape, 
and is temporally advanced by the amount 



r = dd/dw. 



(i) 



If the slope is positive, the wavepacket peak transmits in advance of 
its arrival, whereas if the slope is negative, it transmits with a delay. 

As the wavepacket experiences a sudden temporal shift T g , it also 
experiences a spatial jump, the position of the exiting wavepacket 
being shifted by the amount r\ = c 0 T g , where c 0 is the phase velocity in 
the surrounding medium. Here we present observations of the pro- 
cess of wavepackets passing through a meta-atom with a spatiotem- 
poral jump at —500 Hz by constructing the acoustic meta-atom 
shown in Fig. 2a. Two acoustic resonators A and B are coupled by 
a cavity C, a system designed to emulate an atomic system exhibiting 
electromagnetically induced transparency (EIT) 28,29 . 

Resonator A is made from an 80 mm section of 15.5 mm inner- 
diameter tube with a thin (—10 fim thick) tensioned polyethylene 
membrane at one end. This end of the tube is connected to the 
coupling cavity C of length 51.8 mm and inner diameter 88 mm, 



whereas the other end is open to a one dimensional acoustic 
medium-a duct of diameter 30 mm-so that it interacts with external 
acoustic waves. The quality factor for resonator A, —45, is compar- 
able to that of typical Helmholtz resonators 30 . Resonator B consists of 
a cavity of tunable length —200 ± 50 mm and diameter 32.7 mm, 
and a mass-loaded membrane of the same type as in A. Tunable 
metasurfaces mimicking EIT-like effects have only previously been 
demonstrated for electromagnetic systems 31 " 33 . The Q-factor of the 
resonator B is —140, which is significantly better than that of typical 
acoustic resonators such as resonating ducts and Helmholtz 
resonators 30 . 

We shall derive the dynamics of this system with a lumped- 
element approximation 34 , appropriate here since the characteristic 
geometrical parameters are small compared with our acoustic wave- 
lengths (2—0.68 m). The air in tube A can be considered to move 
uniformly, and its displacement £ A (f) in the direction of the open end 
is given by 



m A t a = ~ k ao £ a - b A I a - S A p c + p ( t) S A , 



(2) 



where p(t) andp c are the pressures at the tube-A mouth and in the 
cavity C, respectively, and the effective mass m A is the sum of the 
masses p 0 S A l' A of the air in the tube and the membrane A. Here, p 0 , S A 
and l A are the density of air, the cross section and the effective length 
of the tube, respectively. (See the Supplementary Information for 
further experimental details.) Also, k A0 and b A are the spring con- 
stant of the membrane and the drag coefficient of the air in the tube, 
respectively. Resonator A emits acoustic waves to the duct, and 
so represents a "radiative" resonator 28 . For convenience we shall 
include the radiation damping into the drag coefficient b A . 
Resonator B has spring constant A; memB 5 , and thus its net spring 
constant is k B0 = /c memB + k ail . The contribution of the elasticity 
of the air in the tube B is /c a [ r — Po^q/ Vb? where V B is the volume of 
this tube. The displacement £ B (t) is given by 



m B ^ B = -k B0 ^ B -htB + S B pc, 



(3) 



where the effective mass m B is the sum of the loaded mass, the 
membrane and the air moving with the membrane, b B is the drag 
coefficient, and S B is the cross section of the tube. Here m A and m B are 
estimated to be 0.0233 and 0.56 g, respectively. The loaded mass on 
the resonator B membrane is estimated to be 0.34 g. Resonator B is 
not directly coupled to the duct, and so represents a "dark" res- 
onator 28 . The pressure p c is related to £ A and £ B as follows: 

Sa£a — $b£b 



pc = B 0 - 



V c 



(4) 



where B 0 is the bulk modulus of air. 

Equations (2)-(4) can be rewritten in the form 

M A l A = -k A <; A + (k c S B /S A )Z B -y A t A +p(t)/S A , 
M B l B = -k B Z B + {k c S A /S B )Z A -y B t B , 



(5) 



where M A = m A /S A , M B = m B /S 2 B are acoustic inertances, 
k A = k A0 /S A + B 0 /V c , k B = k B0 /S 2 B + B 0 /V c , k c = B Q IV c we acoustic 
stiffnesses, y A = b A /S A , y B = b B / S B are acoustic drag coefficients of 
the resonators A and B, and V c is the cavity volume C. These equa- 
tions have the classical form of an EIT system 28 . If, in the uncoupled 
state, A and B are set to have the same resonance frequency, the 
coupling through the cavity C removes the degeneracy and results 
in a split resonance response and a consequent sharp phase disper- 
sion. In practice A and B are set to have similar uncoupled resonant 
frequencies. 

The pressure should be continuous, so p { + p r = p t , where p b p r 
and p t are the pressures of the incident, reflected and transmitted 
waves, respectively. Using the reflection and transmission coeffi- 
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Figure 2 | (a) Experimental acoustic meta-atom consisting of two subwavelength- sized membrane-based resonators A and B coupled via cavity C. 
(b) Experimental setup for transmission-coefficient and pulse-delay measurements. The acoustic waves from the source propagate along the duct, 
through the meta-atom, and are absorbed at the other end. 



cients R = p r /pi and T = pjpb one obtains the relation 1 + R = T. 
Also, the particle -velocity boundary condition is determined by 
equating the volume velocities into and out of the meta-atom: 

Sduct Ui = Sduct U r = Sa %A + Sduct Ut , ( 6 ) 

where S duct is the cross -sectional area of the duct, and the particle 
velocities of the incident, reflected, and transmitted waves are 
defined as u b u r and u t , respectively. From the acoustic impedance 
relations, u { = pi/Z 0 , u r = —p r /Z 0 and u t = p t /Z 0 , Eq. (6) becomes 

T T Sduct Pt 

where Z 0 = ^/lp 0 B 0 ) is the impedance of air. These relations yield the 
acoustic-amplitude transmission coefficient: 

T= 1 — r (8) 

i _i_ Zp Sa Qa 
^ 2 S duct p t 

An analytical expression for T can be obtained from the particle 
velocity \ A by solving the coupled equations of motion for the sys- 
tem, Eqs. (5), under sinusoidal excitation conditions (see the 
Supplementary Information). 

Results 

Equation (8) predicts a phase discontinuity across the meta-atom. To 
investigate this experimentally we measured T(co) = \ T((o)\exp(—i3 ((D)) 
for monochromatic waves, where | T(w)\ and 3(co) are the magnitude 
and phase advance, respectively, over the frequency range 400 to 
600 Hz using the setup shown in Fig. 2b, with microphones placed 
160 mm away either side of the meta-atom. The data for l B = 
195 mm (V B = 163.8 cm 3 ) are shown in Fig. 3a by the dots. The 
solid (red) theoretical curves are obtained from the solutions of Eq. 
(8). These curves, based on measured and fitted parameters (see 
Supplementary Information), agree with the experiment with only 
minor differences, the differences presumably stemming from the 
approximate nature of the lumped- element approximation and the 
effective lengths. 

It is evident that our system can produce phase discontinuities 
— ±1 radian for a single meta-atom, analogous to those seen (up 
to ±7i rad 1117,21 ) in electromagnetic metasurfaces or atomic systems. 
Our system provides a group-time shift, which is important in broad- 
band lensing applications 35 . However, for full control of the waves, 
additional measures may be needed because the span of the phase 
shift 3 in our system is not optimal for redirecting beams. Group 
shifts were measured by comparing the time differences of the peaks 
for the incident and transmitted pulses. Figure 3b shows the results of 
experiments with Gaussian wavepackets of spectral width 3 Hz for 



respectively, the group shift as a function of frequency and the 
pressure amplitude as a function of time at the central frequencies 
506 and 512 Hz, corresponding respectively to the the green solid 
circles p and q in Fig. 3a. The incident and transmitted wavepackets 
are shown in blue (larger pulses) and red (smaller pulses), respect- 
ively. At 506 Hz, the peak of the transmitted pulse is delayed by 
54.0 ms (having the same shape but with reduced height), whereas 
at 512 Hz it is advanced by 40.3 ms. 

To access the full spatiotemporal field, we measured the time 
evolution of the pressure distribution in the tube near the meta-atom 
before and after the wavepacket passes through it. The results are 
shown in the normalized space-time diagrams of Fig. 3c. A discon- 
tinuity is clearly seen at the position z = 0 of the meta-atom. The 
trajectories for incident, reflected and transmitted pulses far (>1) 
from the meta-atom were not measured, but should be straight lines 
with slopes ±c 0 , where c 0 = 340 m/s is the speed of sound. At the 
position of the meta-atom, the space-time trajectory makes temporal 
steps up to —54 ms and 36 ms. These steps occur over a distance 
—0.02 m (—0.03/1), governed by the diameter (15.5 mm) ofthe radi- 
ative resonator. This spatiotemporal observation of a positive 
temporal step provides a first detailed view of the "supersonic" trans- 
mission of pulses through a sub wavelength- thickness meta-atom, 
analogous to that of "superluminal" propagation of light through a 
microscopic interface. Negative group delay here does not imply 
causality violation. It just means that the peak of the transmitted 
wavepacket appears before the arrival of the peak of the incident 
one. As demonstrated in earlier works 36 , this is because the peak of 
the transmitted wavepacket is not causally connected to that of the 
incident one: more energy of the wavepacket is absorbed from the 
rear part than from the front part. 

Group shifts in metasurfaces consisting of resonating meta-atoms 
result from rapid changes in phase shift 3 near the resonance fre- 
quencies. Clearly, for maximum group shift the meta-atoms should 
be designed for maximum |d^/da)|. Our structure exhibits a max- 
imum group shift \z g \ —54 ms, corresponding to the distance |cqT^| 
— 18 m, or 27 L Our work greatly improves on the magnitude of the 
discontinuity compared to the value of 4/1 reported by the time- 
domain measurements of Dolling et al. 26 . 

For the case of optical pulses used for digital communications, the 
pulse widths are approximately equal to the coherence length. When 
metasurface lenses are used to refocus such broadband pulses, the 
path lengths for rays passing through different parts of the metasur- 
face are different. To compensate the path differences it is necessary 
to provide appropriate group shifts x g which vary with position over a 
wide tuning range. We tuned the group shift of our meta-atom for a 
fixed central frequency of 500 Hz by adjusting the resonance fre- 
quency of the dark oscillator, i.e., resonator B, via the cavity volume 
V B , as shown by the results of Fig. 4a. The measured incident (larger, 



SCIENTIFIC REPORTS | 4:4634 | DOI: 1 0.1 038/srep04634 



3 



CD 

'E 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 





P / 




• / 


—Theory 




• Exp. 


V q 




400 450 500 550 
Frequency [Hz] 



600 



400 450 500 550 
Frequency [Hz] 



600 




-0.1 0.0 0.1 
Distance [m] 

0.0 0.5 1 .0 

Normalized pressure amplitude 
[arb. units] 



-0.1 0.0 0.1 
Distance [m] 

0.0 0.5 1 .0 

Normalized pressure amplitude 
[arb. units] 



Figure 3 | (a) Dots: measured magnitude I Tl (co) and phase advance 3 (co) of the transmission coefficient for fixed l B = 195 mm (V B = 163.8 cm 3 ). 
Solid red curves: theoretical fits, (b) Time-domain measurements of the incident (larger, blue pulses) and transmitted (smaller, red pulses) acoustic 
pressure amplitude for the same V B . The Gaussian wavepackets have central frequencies 506 and 512 Hz, corresponding to maximum-negative and 
maximum-positive group shifts. The operating frequencies are marked as the green solid circles p and q in (a), (c) Measured space-time diagrams of the 
normalized pressure amplitude for positive and negative group delay. The volume V B was fixed to 170.5 cm 3 , and a single temporally- Gaussian pulse of 
spectral width 3 Hz was generated and recorded by a microphone scanned along the tube. 



blue) and transmitted (smaller, red) wavepackets are plotted in 
Fig. 4b, and are marked by the green solid squares in Fig. 4a. Here 
the group shift is tunable from —55 ms (slightly better than our 
previous — 54 ms) to 29 ms, corresponding to a path difference of 
—42/1. This range is ample for applications to metasurfaces function- 
ing as thin sonic lenses. 3D plots of | T\ and 3 as a function of fre- 
quency and V B are given in Fig. 4c. The curves for a fixed frequency of 
500 Hz are redrawn in the insets. 

Discussion 

As previously mentioned, negative group delay in electromagnetism 
or acoustics ushers in the possibility of "superluminal" or "super- 
sonic" propagation. For electromagnetic pulses this is experimentally 
well established 37 40 . Researchers sometimes account for such effects 
in metasurfaces from the effective bulk constitutive parameters, by 



regarding single-layer metasurfaces as slabs of metamaterials with 
defined thicknesses. However, problems arise from the artificial 
character of the definition of such thicknesses that obscure, at a 
fundamental level, the localized effects at work at the sample bound- 
aries 10 . Instead of bulk parameters such as permittivity, permeability, 
density or modulus, all that needs to be known to characterize a 
metasurface are the delay and attenuation engendered: in Ref. 10 it 
was shown that metasurfaces are well described by modelling them as 
sheets of zero thicknesses that abruptly change the phases and ampli- 
tudes of the reflected and transmitted waves owing to the scattering 
characteristics of the resonant meta-atoms. Our spatiotemporal 
pressure measurements clearly support this picture. It would be 
interesting to extend this research to acoustic-field spatiotemporal 
imaging of waves passing through metasurfaces, and to electromag- 
netic-field spatiotemporal imaging in analogous EIT-like systems; 
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Figure 4 | (a) Measured group shift at the central frequency 500 Hz as a function of the volume of the resonator B, showing a continuous tuning 
range from —55 to 29 ms, equivalent to a range of c 0 z g = 422. (b) Time-domain measurements of the incident (larger, blue pulses) and transmitted 
(smaller, red pulses) acoustic pressure amplitude at chosen values of V B marked by the green solid squares in (a), (c) Measured amplitude and phase 
advance measured across the meta-atom at various frequencies and volumes V B . The curves for a fixed frequency of 500 Hz are redrawn in the insets. 



several reported electromagnetic metamaterials mimicking EIT are 
of planar form 31 " 33,41 ' 42 . Therefore, this electromagnetic-field imaging 
should be feasible using similar structures. 

In conclusion, we have analytically and experimentally demon- 
strated that meta- atoms producing phase discontinuities for mono- 
chromagnetic waves also produce discontinuities in the space-time 
trajectories of wavepackets. Our tunable experimental design, based 
on coupled radiative and dark acoustic oscillators mimicking an EIT 
system, allows us to observe transmitted pulses at —500 Hz that leave 
the meta-atom with group shifts in the range —55 ms to 40 ms, 
corresponding to record values of the normalized time-shift —28 
< c 0 T g /X < 20 for a room-temperature metamaterial system. A 
lumped- element model accurately describes the underlying disper- 
sion. Observation of the entire process of the wavepackets exiting the 
meta-atom both ahead and behind of their arrivals with space-time 
diagrams accurately reveals the discontinuities represented by the 
group time-shift T g = dd/dco. This work opens the way to practical 
focusing devices for acoustic waves based on 2D arrays 22,23 of our 
meta-atom and to novel applications in beam steering and wavefront 
shaping not only in the audio range but also for analogous systems at 
ultrasonic frequencies. 

Methods 

A 10 W loudspeaker driven by a function generator provides sinusoidal or pulsed 
excitation at 1 m distance from the meta-atom. For the anechoic termination we used 
60 m lengths of 30 mm diameter corrugated tubes attached to both ends of the duct, 
each closed with a commercial air filter for automobiles (Hyundai 28113 37101) to 
prevent residual reflections. 

As shown in Fig. 2b, the acoustic meta-atom is positioned in the center of 
the duct. Resonators A and B and cavity C are fashioned from plastic tubes, and 
resonator B is fitted with a plunger to allow its volume to be tuned. The effective 
masses m A and m B of the resonators A and B are estimated to be 0.0233 and 
0.56 g, respectively. The loaded mass on the resonator B membrane is estimated 
to be 0.34 g. 



For the frequency domain measurements, probe microphones are set halfway up 
the duct, as shown in Fig. 2b. For the spatiotemporal measurements these probe 
microphones are scanned along the tube length to produce the space-time plots. 

The envelopes of the Gaussian wavepackets of spectral width 3 Hz (full width at 
half maximum in the power spectrum) were extracted by taking the magnitudes of the 
Hilbert transform of the waveforms, a widely used demodulation scheme for 
amplitude-modulated signals 43 . Then the curves were fitted using a least-squares 
algorithm to find the peaks, and the group shifts were then obtained from the 
differences. 
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